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We extend a quantum kinetic approach to the description of hadronic showers in space, 

•1— I , 

time and momentum space to deep-inelastic ep collisions, with particular reference to ex- 
periments at HERA. We follow the history of hard scattering events back to the initial 
hadronic state and forward to the formation of colour-singlet pre-hadronic clusters and 
their decays into hadrons. The time evolution of the space-like initial-state shower and the 
time-like secondary partons are treated similarly, and cluster formation is treated using a 
spatial criterion motivated by confinement and a non-perturbative model for hadronization. 
We calculate the time evolution of particle distributions in rapidity, transverse and longi- 
tudinal space. We also compare the transverse hadronic energy flow and the distribution of 
observed hadronic masses with experimental data from HERA, finding encouraging results, 
and discuss the background to large-rapidity- gap events. The techniques developed in this 
paper may be applied in the future to more complicated processes such as eA,pp,pA and 
A A collisions. 



1. INTRODUCTION 



The short-distance physics of isolated hard partonic processes in high energy particle 
collisions is nowadays generally well understood within perturbative QCD, either by calcu- 
lating matrix elements with parton final states, or by parton shower evolution based on the 
QCD renormalization group equation. On the other hand, the long-distance dynamics of 
non-perturbative soft processes and of the confinement mechanism in the process of final- 
state parton-hadron conversion is presently not calculable from first principles, and there- 
fore requires phenomenological model building. Nevertheless, over the past two decades, the 
combination of perturbative QCD calculus and realistic modelling of the non-perturbative 
physics has been developed sufficiently to provide an impressively accurate and predictive 
description of a large class of experimental observables connected with large-momentum 
jets i 

However, with the advent of HERA (ep, possibly eA) and the Tevatron (pp), a new 
regime of QCD at high parton density is opening up, with which one is just beginning to 
come to grips. This regime will be further explorable with the future accelerators RHIC 
(^4^4) and LHC (pp,pA, AA). The common novel feature of these machines is the oppor- 
tunity to study the production and evolution of a system of a large number of partons per 
unit phase space AO, = ArAk, which provides a possible source for new phenomena such 
as non-trivial statistical particle correlations, coherence and interference effects, dissipation 
and collective excitations. Examples of the experimental manifestation of such phenomena 
are: in ep (eA) collisions, an enhanced growth of the parton distributions at small Bjorken 
x @, as well as the observation of diffractive events with large rapidity gaps between target 
and current fragmentation regions ]5|, in eA, pA, AA collisions, events with multiple 
parton scattering |7|, ||], the QCD Landau- Pomeranchuk-Migdal effect || and jet quench- 
ing [|l0|]; in AA collisions, the possible formation of a high-temperature, deconfined parton 
plasma ]iT[ . 

To quantify what we mean by high parton density, consider the hard interaction of 
a probing particle with a hadron or nucleus via a momentum transfer Q 3> Aqqd- The 
probe can be, e.g., a photon (in deep- inelastic scattering) or a parton (in hadronic or nuclear 
collisions). The hard interaction probes space-time distances r ~ 1/Q, thereby resolving a 
density of partons in the probed hadron (nucleon) or nucleus which may be characterized 
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by the number of quark and gluon quanta with a definite value of rapidity y ~ ln(l/x) in 
the transverse plane, p qg oc R^ 2 dN qg /dy ~ (AxfN(x,Q 2 )), where /at(x,Q 2 ) denotes 
the sum of quark and gluon parton distributions in a nucleon, i?j\r the nucleon radius, and 
A the number of nucleons. One can distinguish three regions [12]: (i) r <C 1 fm, p qg <C RJ^ 
— the short-distance, low-density regime of perturbative QCD, (ii) r ~ 1 fm — the non- 
perturbative QCD domain of the complex mechanism of confinement; and (hi) r< 1 fm, 
Pqg ~ R^ 2 — the high-density regime where a dense of parton system is probed at short 
distances, so that perturbative methods may be applied, within a statistical approach. 

There are two extreme ways to penetrate a system of partons with large density at 
short space-time distances: one way is deep-inelastic ep scattering (A = 1) at high energy 
in the region of very small Bjorken id, For instance, at HERA, the extrapolation of 
experimental data implies 30-40(!) gluons in a proton at x ~ 10~ 4 [^]. The other way is 
through collisions of heavy nuclei, in which one can reach high parton densities at not so 
very high energies or small x, due to the large number of overlapping nucleons (A 3> 1). 
This presumably can be achieved at RHIC (x ~ 10 _1 -10 -2 ), and certainly at the LHC 
(x ~ 1CP 3 -10~ 4 ). In particular, at the LHC both the conditions of small x and large A 
may be combined. It is clear that the theoretical study of high-density QCD requires the 
development of new methods by recruiting techniques from relativistic many-body physics, 
the kinetic theory of transport phenomena, renormalization group at finite density (and 
finite temperature), etc.. 

The purpose of the present paper is to start looking at this physics from a space-time 
point of view, and to study the dynamics of high parton densities in deep-inelastic ep scat- 
tering (DIS) in the kinematical region covered by the HERA experiments ZEUS and HI. 
In the light of the detailed hadronic measurements at HERA, which provide information 
about the underlying parton and hadron dynamics, our emphasis is the study of the inter- 
play between perturbative partonic processes and non-perturbative hadron formation. We 
employ a kinetic space-time approach to parton-shower evolution combined with a statis- 
tical model of parton-hadron conversion |l5|] that allows us to follow the time development 
of the particle system in both momentum space and position space, i.e., in 7-dimensional 
phase space d?rd A k. 

The space-time structure of the production and evolution of partonic colour charges and 
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their conversion into 'white' hadrons is the key problem in the dynamics of complex multi- 
parton systems. In the context of high-density QCD, insight into this problem is especially 
important, because the presence of many partons close by in phase space, generated by 
the particle dynamics itself, necessarily causes the propagation and interaction of quanta 
to become non-local and to be correlated statistically in position space and colour space. 
As the system evolves, these conditions will change with time and will in general depend 
on the local density of particles. This is to be contrasted with the familiar translation- 
invariant evolution of well-separated parton jets in empty space, in which case space-time 
correlations are absent or irrelevant, because the jets evolve undisturbed by each other. An 
interesting example of a deviation from unscathed jet evolution has been conjectured to 
occur in e + e" annihilation into hadrons via W + W~ production [i~6| ], where the jets from 
the two VK's overlap and crosstalk, so that the interplay between space-time dynamics and 
the colour flow of close-by partons may lead to a noticeable shift in the experimental W 
mass determination [|j 

Summarizing the above arguments, our motivation in the following is two-fold: 

(i) First, we would like to provide an alternative and supplementary analysis of standard 
non-diffractive DIS events, in order to estimate the 'background' [] to the less well 
understood diffractive, 'large rapidity gap' events. To the extent that other parton 
shower models [ Iq , 19, 21, 2^, Q are generally describing this background well, our 



approach should give the same answer, because our additional space-time information 
which is not contained in previous investigations should not contradict the well-known 
parton evolution in momentum space. On the other hand, the space-time geography 
of non-diffractive events may shed some light on the dynamics of diffractive events, 
which presumably undergo a different space-time development. 

(ii) Secondly, since our approach is in principle designed to be universally applicable to 
high-energy collisions involving lepton, hadron, or nuclear beams, we also see ep colli- 
sions at HERA as a learning ground for future eA (HERA?), pp, pA, AA (RHIC, LHC) 
experiments, whose theoretical description certainly requires knowledge of space-time 



The term 'background' is not to be understood literally, because the contribution of diffractive events 
with a large rapidity gap at HERA is of the order of 10 %, which is still comparably small, although 
experimentally significant M. 
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evolution in order to resolve the complex multi-particle dynamics over the expected 
long collision time scales. 

The paper is organized as follows. In Section 2 we review specific features of DIS at 
HERA, with the primary aim of establishing our nomenclature and notation for kinematic 
variables. In Section 3 we introduce the general concept of our model for the space-time 
development of the hadronic system, recalling relevant aspects of our framework of quan- 
tum multiparticle kinetics, the treatment of the initial state, the space-time development 
of the space-like and time-like parton showers associated with initial- and final-state radi- 
ation, and our spatial criterion for the formation of hadronic clusters and their subsequent 
decays. Section 4 presents our main results, including the time development of the rapidity 
distribution, inclusive hadronic spectra and transverse energy flow. Particular attention 
is paid to the distribution of the mass Mx of the observed hadronic final state in events 
without large rapidity gaps, which reflects the details of our cluster formation mechanism 
and hadronization procedure. 

2. SPECIFIC FEATURES OF DIS AT HERA 

For the purpose of clarity and to define quantities used subsequently, we briefly review in 
this Section some basic notions and kinematics, focusing on the conditions of the ep collider 
HERA, where an electron beam and a proton beam with four momenta p e ^ p = (E, p 2 , Oj_) e ,p 
and 

E e = 27 GeV , E p = 820 GeV , ^ = 296 GeV (1) 

collide head-on. For comparison, in the centre of mass of electron and proton, the energies 
are E e ~ E p = 148 GeV, corresponding to a global shift of the proton rapidity as compared 
to (|l|) from \y p \ = 7.45 to \y p \ = 4.45. 

2.1 Event types 

The physics at HERA may be separated in two classes of event types, illustrated in Fig. 1, 
whose definitions are: 

(i) Non-diffractive events (Fig. la): Here the exchanged virtual photon |^| destroys the 
coherence of the incoming proton by a hard scattering off a quark inside the proton, 
2 In the kinematic region investigated, contributions from Z° exchange can be neglected. 
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and breaks up the proton into coloured subsystems, which are essentially a jet (led 
by the struck quark) and the proton remnant system (consisting of the partons that 
have not taken part in the hard interaction) . This class of events is well described by 
the standard QCD hard scattering picture. 

(ii) Diffractive events (Fig. lb): This class is characterized by an interaction in which 
the proton either remains intact or receives some small internal excitation to become 
a relatively low-mass system, and in which the virtual photon also fragments into 
a relatively low-mass system of particles. This leads generally to experimentally- 
observable large rapidity gaps f\ between the outgoing proton and the rest of the 
produced hadronic system, which may be interpreted as the exchange of a colourless 
object (the 'pomeron') between the photon and the proton. 

For the remainder of this paper we consider exclusively the non-diffractive event type, 
which is describable from first principles in terms of the perturbative QCD parton picture, 
and for which our space-time approach in terms of photon-quark hard scattering, parton 
shower evolution and parton-hadron conversion is applicable as an extension of our previous 
work on e + e _ -collisions O]. The diffractive event type will not be addressed here, since 
it requires specific model extensions which we want to avoid at this point. 

2.2 Kinematics 

The pecularities of the kinematics of DIS in general, and of the HERA facility in particular, 
require a clear specification of which Lorentz frame is chosen - an issue which is especially 
important when dealing with the space-time dynamics. The HERA laboratory frame (= ep 
lab) is the actual experimental setup (c.f. Fig. 2a), in which electron and proton beams 
collide head on, but with beam momenta that differ by more than an order of magnitude. 
This is different from the ep centre- of-mass frame (= ep cms) in which electron and proton 
have equal but opposite momentum, and which is shifted in rapidity as compared to the 
laboratory frame. Most convenient for theoretical analyses, however, is the 7p centre-of- 
mass frame (= jp cms), in which the virtual photon and proton collide head on (c.f. Fig. 
2b). 

Our convention in the following is that frame-dependent quantities generally refer to 
3 The terms 'diffractive events' and 'large rapidity-gap events' are often used synonymously. 
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the ep lab or the ep cms (which, as mentioned after ([!]), is related to the former by a 
trivial shift of the proton rapidity by 2.3 units), whereas Lorentz non-invariant quantities 
which refer to the jp cms are marked by an asterisk. For instance, E and k± represent a 
particle's energy and the momentum transverse to the electron-proton axis in the ep system, 
respectively, while E* denotes the energy in the jp cms and the momentum transverse 
to the photon-proton axis. In either frame, we define the negative z-axis by the proton 
direction. 

Let p e (P) denote the electron (proton) incoming momenta and q the space-like photon 
4-momentum, and define the standard Lorentz invariants for DIS as 

s = (p e + P) , Q = -q , x = — — , y 



2Pq P-Pe 
in terms of these measured momenta, where s is the total invariant mass squared of the 
ep system (|l]), Q 2 specifies the invariant mass of the photon, and x, y are the usual di- 
mensionless Bjorken variables, commonly termed the 'scaling variable' and the 'inelasticity 
parameter', respectively. ^From these definitions, one finds for Q 2 > M 2 that 

Q 2 ^xys, W 2 ^ (P + qf ^Q 2 —, (3) 

x 

where W is the invariant mass of the hadronic system, which equals the total cm energy 
in the ~/p cms. Table 1 familiarizes the kinematic relations among variables x, Q 2 , y 
and W 2 with some numerical examples. Fig. 3 presents schematically the phase-space 
regime spanned by these variables, and emphasizes the region in the x-Q 2 plane which is 
experimentally investigated at HERA. The region of former fixed-target experiments is also 
indicated, corresponding to y < 0.01 and Q 2 ^ 100 GeV 2 . 

For the purpose of relating the kinematic conditions in the ep lab to the experimental 
observables measured or calculated in the "fp cms, we need the Lorentz transformation of 
the particle 4- vectors p^ and p*. For instance, the four-momenta of the incoming proton 
and photon, and of the incoming and outgoing (struck) quark, respectively, are in the ep 
lab (Fig. 2a) given by 

P = (E p , 0, 0, -E p ) 

Q 2 rz ^ Q 2 \ 



1 = \ [ yE e -^-,-^(l-y)Q 2 , 0,-yE. ^, 
p q = (xE p , 0, 0, xEp) 



where Q 3> M p is assumed. On the other hand, in the preferable jp cms (Fig. 2b), the 
corresponding momenta are 

• P.. pl 

[I, 0, 0, 1) 



Pa 



2yE e E p 


^AyE e E p - 


Q 2 


2yE e E p 




^4yE e E p - 


Q 2 


x P* 




2yE e E p 





1 , 0, 0, -1 

2yE e E p ' 



K ^4yE e E p — Q 2 I 1 W p '°'°' l + Ay Q E e E p )- (5) 

The invariant differential cross-section for non-diffractive events with a hard photon 
scattering is the convolution of the elementary photon-quark cross-section with the quark 
and antiquark densities in the struck proton, 

TT- = E^T" Q 2 ), (6) 

dxay *-r dxdy 

where the index % labels the quark and antiquark flavors, with e, and fi(x, Q 2 ) denoting the 

corresponding electric charges and (anti)quark distributions of the proton. The associated 

elementary cross-sections d&i are given to lowest order (i.e., before any QCD radiation) by 

„2 p 2 r , 4« 2 
'era °j 



dn - 2 ™"" f ' \i + (i-y) 2 ) + %^(i-y) 



dxdy yQ 2 ^ V / Q 2 

(2 - y) cos(</>) + (1 - y) cos(20) \ , (7) 



(2 - y) cos(0) + — 



where p± pr im = P±prim(cos(j),sm(fi) is the intrinsic transverse momentum of the primary, 
initial quark or antiquark due to the Fermi motion of the partons inside the proton. In- 
tuitively, one expects the value of p± pr i m to be of the order of the inverse proton radius, 
and it is in fact determined experimentally in hadronic collisions as well as in DIS to be ~ 



400-450 MeV g|]. 
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3. THE MODEL 



3.1 General concept 



The central element in our approach is the use of QCD transport theory [25] and quan 



turn field kinetics [14] to follow the evolution of a generally mixed multiparticle system 
of partons and hadrons in 7-dimensional phase-space d?rd 3 kdk° . We include both the 
perturbative QCD parton-cascade development |2^, |2{| , and the phenomenological 

parton-hadron conversion model which we have proposed previously in Refs. |17|], in 
which we consider dynamical parton-cluster formation as a local, statistical process that 
depends on the spatial separation and colour of nearest-neighbour partons, followed by the 
decay of clusters into hadrons. In contrast to the commonly-used momentum-space descrip- 
tion, in our approach we trace the microscopic history of the dynamically-evolving particle 
system in space-time and momentum space, so that the correlations of partons in space, 
time, and colour can be taken into account for both the perturbative cascade evolution 
and the non-perturbative hadronization. We emphasize that one strength of this approach 
lies in the possible extension of its applicability to the collision dynamics of complicated 
multi-particle systems, as in eA, pA and AA collisions, for which a causal time evolution 
in position space and momentum space is essential. 

The model contains three main building blocks which generically embody high-energy 
collisions involving leptons, hadrons, or nuclei in colliders (for DIS ep collisions, the model 
components are illustrated in Fig. 4): 

a. the initial state associated with the incoming collision partners (the beam particles), 

in particular the phenomenological construction of the hadron (nucleus) in terms of 
quark and gluon phase-space distributions; 

b. the parton cascade development with mutual- and self-interactions of the system of 

quarks and gluons consisting of both the materialized partons from parton showers, 
and the spectator partons belonging to the remnants of the collided beam particles; 

c. the hadronization of the evolving system in terms of parton coalescence to colour-neutral 

clusters as a local, statistical process that depends on the spatial separation and colour 
of nearest-neighbour partons, followed by the decay of clusters into hadrons according 
to the density of final hadron states. 
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Such a pragmatical division, which assumes complex interference between the different 
physics regimes to be negligible, is possible if the respective dynamical scales are such that 
the short-range hard interaction, with its associated perturbative parton evolution, and 
the non-perturbative mechanism of hadron formation occur on well-separated space-time 
scales. For DIS, this condition of validity requires min(W 2 ,Q 2 ) > L~ 2 S> Aq CD , meaning 
that the characteristic mass scale for the "yp hard scattering and parton shower development 
(W 2 , Q 2 , or a combination of the two) is larger than the inverse 'confinement length scale' 
L c = 0(1 fm) separating perturbative and non-perturbative domains. Specifically, for DIS, 
it is apparent from (jH) that in the small- x regime probed at HERA (1(T 4 < x < 1(T 3 ), one 
has 60 < W <300 GeV for 10 < Q 2 < 300 GeV 2 , so that the above requirement is well 
satisfied. We emphasize however, that in our model the interplay between perturbative 
and non-perturbative regimes is controlled locally by the space-time evolution of the mixed 
parton-hadron system itself, rather than by an arbitrary global division between parton 
and hadron degrees of freedom. 

We now turn to the specific case of DIS, and in the following subsections we will discuss 
the above components in more detail. 

3.2 Framework of quantum kinetics for multiparticle dynamics 

^From quantum kinetic theory, one can obtain a space-time description of multiparticle 
systems in high-energy QCD processes, as has been discussed formally in Ref. [[l4||. Applied 
to the concept of our model, as outlined in Sec. 3.1, this framework allows us to express 
the time evolution of the mixed system of incoherent partons, composite clusters, and 
physical hadrons in terms of a closed set of integro-differential equations for the local 
phase-space densities of the different particle excitations. The definition of these phase- 
space densities ('Wigner densities'), denoted by F a , where a = p,c,h labels the species of 
partons, prehadronic clusters, or hadrons, respectively, is: 



particles of type a at time t with position in r + df, momentum in k + dk, and energy in 
E + dE (or equivalently invariant mass in k 2 + dk 2 ). The F a are the quantum analogues of 
the classical phase-space distributions, and contain the essential microscopic information 




dN a (t) 



(8) 



d?rd?kdE 



where k 2 = E 2 — k 2 



can be off or on mass shell. The densities (H) measure the number of 
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required for a statistical description of the time evolution of a many-particle system in 
complete phase space, thereby providing the basis for calculating macroscopic observables 
in the framework of relativistic kinetic theory. 

The Wigner densities (|J) are determined by the self-consistent solutions of a set of 
transport equations (in space-time) coupled with renormalization-group-type equations (in 
momentum space). Referring to Refs. fl7| , |l5| , 14] for details, we remark that these equations 



can be generically expressed as convolutions of the densities of radiating or interacting 
particles Fp with specific cross sections Ij for the processes j, yielding the following closed 
set of balance equations for the space-time development of the densities of partons F p , 
clusters F c and hadrons Fh, 

k ■ d r F p (r, k) = F p , o J(p' -> pp") - F p o J(p -> p'p") - F p F p , o I(pp' -» c) (9) 
k ■ d r F c (r, k) = F p F p , o I(pp -> c) - F c a J(c -» fc) (10) 
k-d r F h {r,k) = F c oI(c^h), (11) 

where k ■ d r = k^d/dr^. We remark that eq. @ implicitly embodies the momentum 
space (k 2 ) evolution of partons through the renormalization of the phase-space densities 
F p , determined by their change k 2 dF p (r,k)/dk 2 with respect to a variation of the mass 
(virtuality) scale k 2 in the usual QCD evolution framework [26, 27] 0. Each of the terms on 



the right-hand side of corresponds to one of the following categories (c.f. Fig. 4): 

(i) parton multiplication through radiative emission processes on the perturbative level, (ii) 
colourless cluster formation through parton recombination depending on the local colour 
and spatial configuration, (iii) hadron formation through decays of the cluster excitations 
into final-state hadrons. Each convolution F o I of the density of particles F entering 
a particular vertex I includes a sum over contributing subprocesses, and a phase-space 
integration weighted with the associated subprocess probability distribution of the squared 
amplitude. 

The equations (p|)-([Tl|) reflect a probabilistic interpretation of QCD evolution in space- 
time and momentum space in terms of sequentially-ordered interaction processes j, in which 

the rate of change of the particle distributions F a (a = p, c, h) in a phase-space element 
4 For pre-hadronic clusters and hadrons, we assume renormalization effects to be comparatively small, 
so that their mass fluctuations Ak 2 /k 2 can be ignored to first approximation, implying k 2 dF c (r,k)/dk 2 = 
k 2 dF h (r,k)/dk 2 = 0. 
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d?rd k is governed by the balance of gain (+) and loss (— ) terms. The left-hand side 
describes free propagation of a quantum of species a, whereas on the right-hand side the 
interaction kernels I are integral operators that incorporate the effects of the particles' 
self and mutual interactions. This quasi-classical, probabilistic character of high-energy 
particles is essentially an effect of time dilation, because in any frame where the particles 
move close to the speed of light, the associated wave-packets are highly localized to short 
space-time extent, so that long-distance quantum interference effects are generally very 
small. 

3.3 Scheme of solution and choice of Lorentz frame 

In the above kinetic approximation El] to the multi-particle dynamics, the probabilistic 



character of the evolution equations (O)-(ll) allows one to solve for the Wigner densities 



F a (r, k) by simulating the dynamical development as a Markovian process causally in time. 
Because it is an initial-value problem, one must specify some physically appropriate initial 
condition F a (to,r,k) at starting time to, such that all the dynamics prior to this point is 
effectively embodied in this initial form of F a . The set of kinetic equations (p|)-(|l~ll) can 
then be solved in terms of the evolution of the Wigner densities F a for t > to using Monte 
Carlo methods to simulate the time development of the mixed system of partons, clusters, 



and hadrons in position and momentum space [17, 25 1. 

In the next subsections we explain in more detail the different components for the case 
of DIS, namely, the initial-state Ansatz, parton shower development, and parton-hadron 
conversion. The overall concept of the simulation is illustrated in Fig. 4 and can be sum- 
marized as follows: given the initial state of the photon and the proton disassembled into 
its parton content, the hard interaction of the photon with one of the quarks occurs at time 
t = 0. Specifying the initial state at some earlier time to < 0, and with the hard scattering 
variables chosen from the cross-section, the phase-space distribution of particles at t = can 
be calculated and then evolved in small time steps forward, until stable final-state hadrons 
are left as freely-streaming particles. The size of time steps is chosen as At = O(10~ 3 fin), 
so that an optimal resolution of the particle dynamics in space and energy-momentum is 
achieved. The partons propagate along classical trajectories until they interact, i.e., decay 
(branching process) or recombine (cluster formation). Similarly, the clusters so formed 
travel along classical paths until they convert into hadrons (cluster decay). The corre- 
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sponding probabilities and time scales of interactions are sampled stochastically from the 



It is clear that the description of particle evolution is Lorentz-frame dependent, and 
a suitable reference frame must be chosen (not necessarily the laboratory frame). When 
computing Lorentz-invariant quantities, such as cross sections or final-state hadron spectra, 
the particular choice is irrelevant, whereas for non- invariant observables, such as energy 
distibutions or space-time-dependent quantities, one must at the end transform from the 
arbitrarily-chosen frame of theoretical description to the actual frame of measurement. 
Furthermore, at HERA even experimental analyses are often carried out in the 7p cms 
(||), rather than the ep lab (|J). For our purposes it is most convenient to choose the 
overall centre-of-mass frame of the colliding electron and proton, the ep cms, as the global 
frame with respect to which the evolution of the collision system is followed []. Recall our 
convention that the ep collision axis defines the z-axis, with the electron (proton) moving 
in the positive (negative) z direction. The incoming 4-momenta p e and P involve therefore 
no transverse components, and are ^ 



where P cm = W s — (m e + M p ) 2 Js — (m e — M p ) 2 / (2i/i) is the ep cm momentum. 



3.4 Initial state 

The incoming electron is considered as a point-like object carrying the full beam energy, 
meaning that we neglect any QED or QCD substructure of the electron, as well as initial- 
state photon radiation by the electron. We assume that the electron emits the virtual 
photon of invariant mass Q 2 = — q 2 at time t = — Q _1 , so that t = characterizes the point 
when the photon hits the incoming proton, as is depicted in Fig. 4. 

The incoming proton, on the other hand, is decomposed into its parton substructure by 

phenomenological construction of the momentum and spatial distributions of its daughter 

5 However, to make contact with the HERA experiments, most of our results will be discussed later in 

the 7p cms rather than the ep cms, unless specified otherwise. 

6 We emphasize that in our calculations we use exact kinematics, and take into account proton, electron 

and quark masses. 



relevant probability distributions in the kernels I of eq. (|9|)-(|Tl|). 



Pe 



P 




(12) 
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partons on the basis of the experimentally-measured proton structure functions and elastic 
proton form factor. Here it is important to distinguish between the scales Q 2 and Qq 
(c.f. Fig. 4): The hard scattering scale Q 2 = —q 2 is set by the momentum transfer q 
between electron and proton and determines the parton structure as seen by the virtual 7 
after the initial state radiation of the struck quark. The initial resolution scale Qq, on the 
other hand, determines how detailed the parton phase-space density in the proton would be 
resolved before the initial state radiation. Hence, in accord with (||), we introduce the initial 
parton phase-space distribution (r, p) as the number density of partons in a phase-space 
element d?rd?pdE at time t = to within in the proton at an initial resolution scale Qq = 1 
GeV 2 . We assume the following factorized form: 



F®(r,p) = F a (r,p) 



t=t , | p 2i~r>2 



P a (p,P;Q 2 )oR a (p,f,R) 



i) 



(13) 

t=t 



The right-hand side P a o R a is a convolution of an initial momentum distribution P a and a 
spatial distribution R a , with the subscript a = g,q%,qi labeling the parton species (gluons 
or (anti)quarks of flavor i = 1, . . . n/). The 4-vectors p = p^ = (E,p) and r = r^ = (t,r) 
refer to the partons, whereas P = (0,0, —P C m) and R = refer to the initial 3-momentum 
and the position of the parent proton at t = to in the ep cms. The partons' energies 



E = \Jp 2 — Qq take into account initial space-like virtualities p 2 < 0, which reflects the 
fact that before the collision the partons are confined inside the parent proton and cannot 
be treated as free particles (meaning that they do not have enough energy to be on mass 
shell, but are space-like off-shell). The initial momentum distribution is taken as 

Pa(p,P;Ql) = (I) f a (x,Ql) g(p±) 5 (p z - ^ 5 2 (P ± ) . (14) 
Here x and x are the partons' longitudinal momentum and energy fractions, respectively, 



Pz Ep 



f +i ^r (15) 



and the functions f a (x, Qq) are the usual (measured) quark and gluon structure functions of 

the proton which specify the longitudinal momentum distribution, whereas the transverse 

momentum distribution g(p±) = (27rpo) _1 exp [— p\/pq\ takes into account the uncertainty 
7 We use the GRV structure function parametrization which describes quite accurately the HERA 
data even at low Q 2 and very small x. 
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of the transverse momentum ("Fermi motion") due to the fact that the initial partons are 
confined within the proton. The latter is inferred from experimental analyses |24f| , with 
Po = 0.42 GeV, corresponding to the mean primordial transverse parton momentum 
The normalization is such that 

rl poo 

E / dxxf a (x,Q 2 ) = 1 / d 2 p ± g(p ± ) = 1 (16) 

EP a (p,P;Q 2 ) = n{P,Ql) , (17) 



(2nf(2E) 

where n(P, Qq) has dimension 1/volume and gives the total number density of partons in 
the proton with momentum P, when resolved at the scale Qq. Finally, we impose the 
constraint that the total invariant mass of the partons equals the proton mass M p , 

X>J - feft*) - fix) - (x>.) = M l > (is) 

where the summation j = 1, . . . n(P, Q 2 ) runs over all partons resolved at Q 2 , as constrained 
by (|l6|) and (|l7]). With the partons' 3- momenta determined from the distributions in x and 
p±, the requirement (|l8|) fixes the relation between energy and momentum by assigning 
to each parton an initial space- like virtuality such that p 2 = E 2 — p 2 < 0. With this 
prescription, the resulting distribution in p 2 is approximately Gaussian with a mean value 



of \J (p 2 ) ~ 500 MeV, i.e. the typical initial virtuality of the partons is about Qo/2. 
3.5 Parton cascade development 

With the above construction of the initial state in terms of the incoming electron and 
photon, and the parton cloud of the proton, the dynamical development of the system can 
now be traced according to the kinetic equations (]9|)-(|il"|), starting from t = 0. In our 
statistical picture, the initial-state parton ensemble represents a particular fluctuation of 
the proton wave function that has developed between t = to ~ Qq 1 < and time t = 0, 
at which the photon with resolution Q 2 picks according to the cross-section @ a quark 
with specific flavour and momentum p = xP out of the incoming parton cloud, while the 
other partons are viewed as unaffected by the short-range "fq interaction. Consequently, as 
illustrated in Fig. 4, the early stage of the time evolution is characterized by two different 
physics elements: a) the parton showers initiated by the quark that is struck out of the 
original proton wave function through the momentum transfer from the virtual photon, 
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and b) the propagation of the remnant system consisting of the other initial partons, that 
remain spectators of the hard process and form the coherent remnant of the original proton. 



For the parton shower development we employ the well-established jet calculus [28, ^] 
based on the 'modified leading logarithmic approximation' (MLLA) to the QCD evolution 
of hard processes [M, E7| . A parton shower then reduces to a strictly-ordered sequence of 
elementary branchings q — > qg, g — > gg, g — > qq, which can be described stochastically as a 
Markov cascade in position and momentum space. We distinguish initial-state, space-like 
branchings of the selected quark before it reaches the "fq vertex, and final-state, time-like 



radiation off the struck quark after the hard r yq interaction [23, 31] q. The separation 
into two 'hemispheres' divided by the ^q vertex is illustated in Fig. 5: it refers to both 
the chronological order along the real time axis and to the order of emission vertices in 
momentum space. The initial quark that is picked out by the photon evolves from the 
remote past t = to < towards the hard interaction by sequential branchings pj — ► Pj+i + 
p'j + i, in each of which one of the daughters continues with increasing space- like virtuality 
\Pj+i\ > \Pj\ (where p 2 ,p 2 + i < 0), while the other one acquires a time-like virtuality p' 2 + \ > 
and may develop a time-like shower of its own. The space-like shower is consequently 
characterized by increasing virtualities \p 2 \ < , decreasing energies and increasing 

opening angles, as the quark approaches the hard vertex at t = with \p^\ ~ Q 2 - Once 
the evolved quark has been struck by the photon, the momentum transfer provides the 
outgoing quark with enough energy-momentum to become a real excitation at t = and to 
obtain a time-like virtuality ~ Q 2 . This materialized quark initiates now a shower of 
sequential time-like branchings k m — > k m ~\ + k' m _± in which both daughters are time-like 
(i.e. ^-l'^TO-i > 0) w ith decreasing virtualities fc^-i < ^m> decreasing energies and 
decreasing opening angles. The branching chain continues into the remote future until it is 
terminated by the hadronization, which we model as the coalescence of neighboring partons 
in a cascade, followed by conversion to hadrons (Sec. 3.7 below). 

The specific feature of our approach is that, in addition to the definite virtuality and 
momentum, each elementary vertex has a certain space and time position which is obtained 

by assuming that the partons in the shower propagate on straight-line trajectories in be- 
8 This separation implies the neglect of interference between the initial- and final-state showers, a common 
conceptual defect that is approximately cured by matching on to the lowest-order 0(a s ) matrix element. 
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tween the branchings. In the MLLA framework, the basic properties of both space-like and 



time-like showers are determined by the DGLAP equations [32|, but with essential differ- 
ences in time ordering, kinematics and the treatment of infrared singularities associated 
with soft gluon emission. 

3.5.1 Space-like parton shower 

As mentioned above and depicted in Fig. 5, the space-like cascade starts at some time 
t = to ~ — Qq 1 before the actual hard scattering at t = 0, with the initiating parton of 
virtuality \pq\ ~ Qq = M 2 embodied in the parton cloud of the incoming proton, and 
proceeds up to p 2 = p 2 n ~ — Q 2 at the hard vertex set by the space-like photon virtuality. 
The emitted partons on the side branches, on the other hand, are not connected directly 
with the 7<7 vertex, but evolve independently as time-like quanta. In the cascade sequence 
both collinear and soft coherent branchings are properly included [j33|, if the development of 
the chain is described in terms of 'angular-ordering' variables (rather than the virtualities 

P% 

ft) '■:]<,■: • Cj+i = P °'5 +1 - l-cosfl ,j+i (0<j<n), (19) 

E E j+1 

where pj = (Ej,pj) and p'j + \ = (Ej +1 ,pj +1 ) are assigned as in Fig. 5 for the j th branching 
Pj ~^ Pj+iPj+i- The space-like cascade is then strictly ordered in the variable p 2 + \ > ft 2 , 
which is equivalent to the ordering of emission angles, Ej6 0: < Ej + i9 0t ■ 

Because the presence of the external hard interaction at t = and Q 2 sets a physical 
boundary condition on the kinematical evolution of the cascade, it is technically advanta- 
geous to reconstruct the cascade backwards in time starting from t = at the hard vertex 
Q 2 and trace the history of the struck quark back to Qq at t = to . The method used here 



is a space-time generalization of the 'backward evolution scheme' [31, 34]. To sketch the 
procedure, consider the space-like branching p n _i — > p n p' n which is closest to the jq vertex 
in Fig. 5. The virtualities satisfy M] \p%\ > \Pn-i\i an d Pn^Pn-l < (space-like) but 
p' 2 > (time- like). The relative probability for this branching to occur between ft 2 and 
ft 2 + dft 2 is given by 

x { F[r :- 1]Xn ^f ) \ T^At), (20) 
^ F(r n ;x n ,ft 2 ) J K n K ' 
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where Xj = (pj) z /P z (j = n, n — 1) are the fractions of longitudinal proton momentum P z , 
with F(rj) Xj,p 2 ) = F(rj,pj) the corresponding parton distriibutions introduced before, 
and the variables 

En x n E„ x n —i x n , , 

1 — z = — — ~ (21) 



Efi—X %n—l E n —i X n —i 

specify the fractional energy or longitudinal momentum of parton n and n' , respectively, 
taken away from n — 1. The function a s /(2irp 2 ) 7(2) is the usual DGLAP branching prob- 
ability in the MLLA, with 7(2;) giving the energy distribution in the variable z. The last 
factor in (|20| ) determines the time interval in the ep cms, At = t n — i n _i, that is associated 
with the branching process n — 1 — > nn' . We take here simply 

T^(At) = S ^'^ P, - At) , (22) 

which accounts for the formation time of n by its mother re-Ion the basis of the uncertainty 
principle: At = AE/\pl\, AE ~ (x n - x n _i) P z . 

The "backwards evolution" of the space-like branching p n -i — ► Pn+ p' n is expressed in 
terms of the probability that parton (n — 1) did not branch between the lower bound 
given by the initial resolution scale Qq, and p 2 . In that case, parton n can not originate 
from this branching, but must have been produced otherwise or already been present in 
the initial parton distributions. This non-branching probability is given by the Sudakov 
form-factor for space-like branchings: 

S n (x n ,p 2 ,p 2 ; At) = e X pj-]r [ P [ Z+iP) dV^ix^z^p' 2 ; At)) , (23) 

I a J P J z-(P') J 

where the sum runs over the possible species a = g, q, q of parton n—1. The upper limit of 
the p 2 -integration is set by p 2 ^Q 2 , associated with the scattering vertex of quark n with 



the photon in Fig. 5. The limits z± are determined by kinematics [33]: Z-(p) = Qo/p 
and z+(p) = 1 — Qo/p. The knowledge of S n (x n ,p 2 , Qq) is enough to trace the evolution 
of the branching closest to the hard vertex backwards from p\ at t = t n = to Pn_i 
at i n _i = — Xn/lPnl P z . The next preceding branchings p n -2 — > p n -ip' n -i, etc., are then 
reconstructed in exactly the same manner with the replacements t n — > t n —\, x n — > x n —\, 
Pn ~* Pn-i-> an d so forth, until the initial point p^ at to = —Qq 1 is reached. 
3.5.2 Time-like parton shower 
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Time-like cascades are initiated by the secondary partons, i.e., those that emerge from 
the side branches of the initial-state radiation from the scattering quark before t = 0, 
as well as those that are produced by final-state emission from the scattered quark after 
t = 0. Consider the time- like cascade in Fig. 5, that is initiated by the outgoing quark of 
momentum k = k m emerging from the hard vertex at t = and off shell by an amount 

Again an angular-ordered (rather than virtuality-ordererd) time evolution of the cascade 
is employed to incorporate interference effects of soft gluons emitted along the tree in Fig. 
5. In contrast to ([TJ]), the time-like version of the angular evolution variable is [35| 



k) = E] , = P ' 1 ' V J 1 ~ 1 - cos 9 {rlUriy (m > j > 1) . (24) 

'\i i j i 

so that the time-like cascade can be described by a fc 2 -ordered (rather than /c 2 -ordered) 
evolution, which corresponds to an angular ordering with decreasing emission angles Ojji > 

Proceeding analogously to the space-like case (c.f. (§g)), the probability dV^> m _ x for 
the first branching after the jq vertex, k m — ► k m -\k' m _i with k^-i, is given by the 



space-time extension |14j, |2JJ] of the usual DGLAP probability distribution |32j, 



dV^m-i^; At) = ^ Z ^7MHM'W T^(At) , (25) 

where T^ T \At) is the probabilty that parton m with virtuality and corresponding 
proper lifetime r m oc decays within a time interval At, 

im ^ = 1 --(-^)- (26 » 

The actual lifetime of the decaying parton m in the ep cms is then t m (k) = j/r m (k), where 
t q (k) ~ 3S/(2a s A; 2 ) for quarks and t g (k) « E/(2a s k 2 ) for gluons (8|. As before, i*j denotes 
the local density of parton species j = rn,m — l, and a s / (2^)7(2) is the DGLAP branching 
kernel with energy distribution 7(2). The probability (^) is formulated in terms of the 
energy fractions carried by the daughter partons, 

z = %^ , 1-* = %± , (27) 
with the virtuality fc m of the quark m related to z and £ through k 2 m = k 2 ni +k' r %_ 1 +2E 2 a z(l— 



z)£, and the argument k in the running coupling a s in (25) is 33] k = 2z (1 — z) E^ 
k\. 
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The branching probability fl2q ) determines the distribution of emitted partons in both 
coordinate and momentum space, because the knowledge of four-momentum and lifetime 
(or At between successive branchings) give the spatial positions of the partons, if they are 
assumed to propagate on straight paths between the vertices. The probability that parton 
m does not branch between k and a minimum value = /Xq is given by the exponentiation 
of (|2q), yielding the Sudakov form-factor for time-like branchings: 

rk 2 ^ rz + (k') 

which is summed over the species a = g,q,q of parton m — 1. The integration limits fcg 
and z± are determined by the requirement that the branching must terminate when the 
partons enter the non-perturbative regime and begin to hadronize. As we discuss later, 
this condition can be parametrized by the confinement length scale L c = 0(1 fm) with 
k~Q 2 > L~ 2 = //q, and z + (k m ) = 1 — z_(k m ) = [io/\J ^k m , so that for z+Qzq) = 2_(fco) = 1/2 
the phase space for the branching vanishes. 



T m (k 2 ,k 2 ;At) = exp{- C E f Z+ik) d^m-x^k' 2 ; At)) , (28) 



^0 a 



The Sudakov form factor (28) determines the four- momenta and positions of the partons 
of a particular emission vertex as we sketched above for the first branching, but subsequent 
branchings are described completely analogously by replacing t m — > t m _i, x m — ► x m -%, 
km ~^ k m -i, etc.. Hence T(k 2 , kq] At) generates the time-like cascade as sequential branch- 
ings starting from t = at the hard vertex forward in time, until the partons eventually 
hadronize as discussed below. 

3.6 Cluster formation and hadronization 

Both the cluster formation from the collection of quarks and gluons at the end of 
the perturbative phase and the subsequent cluster decay into final hadrons consist of two 
components: 

(i) The recombination of the secondary time-like partons, their conversion into colourless 

parton clusters and the subsequent decay into secondary hadrons. 

(ii) The recombination of the primary space-like partons that remained spectators through- 

out the collision development into beam clusters and the fragmentation of these clus- 
ters. 
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The important assumption here is that the process of hadron formation depends only 
on the local space-, time-, and colour-structure of the parton system, so that the hadroniza- 
tion mechanism can be modelled as the formation of colour-singlet clusters of partons as 
independent entities (pre-hadrons), which subsequently decay into hadrons. This concept 
is reminiscent of the 'pre-confinement' property j3(| of parton evolution, which is the ten- 
dency of the produced partons to arrange themselves in colour-singlet clusters with limited 
extension in both position and momentum space, so that it is suggestive to suppose that 
these clusters are the basic units out of which hadrons form. 
3.6.1 Cluster formation 

(i) Parton clusters: Parton clusters are formed from secondary partons, i.e. those 
that have been produced by the hard interaction and the parton shower development. The 
coalescence of these secondary partons to colour-neutral clusters has been discussed in detail 
in Refs. [15, so that we confine ourselves here to the essential points. Throughout the 



dynamically-evolving parton shower development, we consider every parton and its nearest 
spatial neighbour as a potential candidate for a 2-parton cluster, which, if colour neutral, 
plays the role of a 'pre-confined' excitation in the process of hadronization. Within each 
single time step, the probability for parton-cluster conversion is determined for each nearest- 
neighbor pair by the requirement that the total colour charge of the two partons must give 
a composite colour-singlet state, and the condition that their relative spatial distance L 
exceeds the critical confinement length scale L c . We define L as the Lorentz- invariant 
distance Ly between parton i and its nearest neighbor j: 

L(n,rj) = = min(Aji, . . . , Ay, . . . , A in ) , (29) 



where Ay = yVyry )M , r%j = — rj, and the probability for the coalescence of the two 
partons i, j to form a cluster is modelled by a distribution of the form 

Il^-c oc (l — exp (-AF Lij)\ ~ 1 - exp ( % - ) if L < Ly < L c , (30) 



L c -L 



if Lq < Lij < L c , and (1) if Ly < Lq (Lij > L c ). Here AF is the local change in the free 
energy of the system that is associated with the conversion of the partons to clusters, and the 
second expression on the right side is our parametrization in terms of Lq = 0.6 fm and L c = 
0.8 fm that define the transition regime. As we studied in Ref. [17], the aforementioned 



20 



colour constraint, that only colourless 2-parton configurations may produce a cluster, can 
be incorporated by allowing coalescence for any pair of colour charges, as determined by the 
space-time separation Lij and the probability (|30|), however, accompanied by the additional 
emission of a gluon or quark that carries away any unbalanced net colour in the case that 
the two coalescing partons are not in a colourless configuration. 

(ii) Beam clusters: The remaining fraction of the longitudinal momentum and energy 
that has not been redirected and harnessed by the interaction with the photon is carried 
by the primary partons of the initial proton, which remained spectators throughout. In 
our approach these partons maintain their originally assigned momenta and their space- 
like virtualities. Representing the beam remnant, they may be pictured as the coherent 
relics of the original proton wavefunction. Therefore the primary virtual partons must be 
treated differently than the secondary partons which are real excitations that contribute 
incoherently to the hadron yield. In the ep cms the primary partons are grouped together 
to form a massive beam cluster with its four-momentum given by the sum of the parton 
momenta and its position given by the 3-vector mean of the partons' positions. 
3.6.2 Hadronization of clusters 

(i) Parton clusters: 

For the decay of each parton cluster into final-state hadrons, we employ the scheme pre- 
sented in Refs. [^, 17, 38]: If a cluster is too light to decay into a pair of hadrons, it is taken 



to represent the lightest single meson that corresponds to its partonic constituents. Other- 
wise, the cluster decays isotropically in its rest frame into a pair of hadrons, either mesons 
or baryons, whose combined quantum numbers correspond to its partonic constituents. The 
corresponding decay probability is chosen to be 

„ ., f mc dm ( m\ , , 

U c ^ h = T c (E c ,m 2 c ) M exp , (31) 



where M is a normalization factor, and the integrand is a Hagedorn spectrum [37| that 
parametrizes quite well the density of accessible hadronic states below m c which are listed 
in the particle data tables, and mo = m n . In analogy to (|26|), T c is a life-time factor giving 
the probability that a cluster of mass m\ decays within a time interval At in the global 
frame, here the ep cms, 

jae ^ = 1 - exp (-aitW) ■ <32) 
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with the Lorentz-boosted life time t c = 7 c t c ~ E c /m 2 . In this scheme, a particular cluster 
decay mode is obtained from (|3l]) by summing over all possible decay channels, weighted 
with the appropriate spin, flavour, and phase-space factors, and then choosing the actual 
decay mode acording to the relative probabilities of the channels. 

(ii) Beam clusters: 

The fragmentation of the beam cluster containing the spectator partons mimics in 
our model what is commonly termed the 'soft underlying event', namely, the emergence 
of those final-state hadrons that are associated with the non-perturbative physics which 
underlies the perturbatively-accessible dynamics of the hard interaction with parton shower 
fragmentation. 

In the spirit of Ref. [33], we employ a (suitably modified for our purposes) version 



of the soft hadron production model of the UA5 collaboration [39|, which is based on a 
parametrization of the CERN pp collider data for minimum-bias hadronic collisions. The 
parameters involved in this model are set to give a good agreement with those data. 



We view soft hadron production as a universal mechanism [40] that is common to all 
high-energy collisions that involve beam hadrons in the initial state, and that depends 
essentially on the total energy- momentum of the fragmenting final-state beam remnant. 
Accordingly, we assume that the fragmentation of the final-state beam cluster depends 
solely on its invariant mass M, and that it produces a charged- particle multiplicity with a 



binomial distribution [39] 



Pin) - T{n + k) {Tl/ky (33) 
{> n\T(k) (l + n/k) n + k ' { ' 



where the mean charged multiplicity n = n(M 2 ) and the parameter k = k(M 2 ) depends on 
the invariant cluster mass [] according to the following particle data parametrization [|39| . 



n(M 2 ) = 10.68 (M 2 ) - 115 - 9.5 k{M l ) = 0.029 ln(M*) - 0.064 . (34) 



Adopting the scheme of Marchesini and Webber |33|| , the fragmentation of a beam cluster 
of mass M proceeds then as follows: First, a particle multiplicity n is chosen from 



9 Notice that in our model M fluctuates statistically, as a result of fluctuations of the initial-state 
parton configuration in the proton and the variation of the hard scattering variables x and Q 2 . Hence the 
distribution (^) and the mean multiplicity ([ji]) vary from event to event. This is in contrast to the original 
UA5 model, in which the fixed beam energy \fs/2 controls the energy dependence of soft hadron production. 
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and the actual charged particle multiplicity is taken to be n plus the modulus of the beam 
cluster charge. Next, the beam cluster is split into sub-clusters (qi^), (Q2Q3), ■ ■ ■ (it = u, d), 
which are subsequently hadronized in the beam cluster rest frame, in the same way as 
the parton clusters described in the preceding subsection. To determine the sub-cluster 
momenta, we assume a mass distribution 

P{M) = c (M — 1) exp [-a(M - 1)] , (35) 

with c a normalization constant and a = 2 GeV -1 , resulting in average value of (M) « 1.5 
GeV. The transverse momenta are taken from the distibution 



P(P±) = c' P± exp 



bJp\ + M 2 



(36) 



with normalization d and slope parameter 6 = 3 GeV -1 , and the rapidities y are drawn 
from a simple flat distribution P(y) oc const, with an extent of 0.6 units and Gaussian 
tails with 1 unit standard deviation at the ends. Finally, all hadronization products of 
the sub-clusters are boosted from the rest frame of the original beam cluster back into the 
global frame, i.e. the ep cms. 



4. MODEL RESULTS FOR NON-DIFFRACTIVE DIS AT HERA 

4.1 Characteristic evolution of small- cc versus large-a: scattering events 

The kinematics of DIS has very different consequences in the small-x and large-x regime, 
as we shall discuss now within our model. Specifically, we distinguish here and in the 
following the two distinct regimes 

'small' x: 1.7 • 1(T 4 < x < 2.3 • 1(T 3 'large' x: x > 5 • 1(T 3 , (37) 

where the small-x regime is the typical range probed at HERA, and part of the large-x 
range (x •> 5 • 10~ 2 ) corresponds to previous fixed-target experiments (c.f. Fig. 3). Table 
1 provides the corresponding mass W of the hadronic system (equal to the total hadronic 
energy in the jp cms), where the small-x regime is the HERA range which we primarily 
focus on in the following. 

Figure 6 illustrates vividly the differences between the small-x (left panel) and large-x 
(right panel) kinematics for typical HERA values of Q 2 = 8/14/28 GeV 2 . The top plots 
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show the associated probability distributions for the occurrence of a particular mass W of 
the hadronic system produced by the hard interaction in the specified Q 2 - and x-range. As 
could be expected, the differences between the two x-domains are striking: not only the 
shape, but also the mean values of the distributions are very distinct (note the different 
scales of the W axis). The most probable W values lie between 100-200 GeV (small x) 
and 10-20 GeV (large x). The plots in the middle show the W dependence of the total 
hadron multiplicity calculated within our model. The shape of the curves is again rather 
different, which is a direct consequence of the probability distribution in W, and the phase 
space available for a given W. On the other hand, the total numbers of produced hadrons 
(Nh ~ 35 — 40) in events around the most probable W are very similar. The bottom plots 
show the corresponding mean values of the xf = 2p z /W, with p z in the ep cms along the 
beam axis in the opposite direction to the incoming proton, i.e., the fractional longitudinal 
momentum carried by the final-state hadrons which emerge from the parton shower and 
fragmentation of the struck quark jet. Again the W dependences of (xp) are very different 
in the small-x and large-x regimes, with the typical xf ranges 0.03 < (xp) ^ < 0.13 and 
(xp) > 0.5, respectively. 

The distinct characters of small-x and large-x DIS events in the kinematic regimes 
discussed above are accompanied by different space-time evolution patterns of the particles 
in position and momentum space. As explained in Sec. 3, our approach allows us to 
follow the time rate of change of the particle densities and associated spectra. Fig. 7 
exhibits the time evolution of the rapidity (y) distribution, and the particle distributions 
in longitudinal (z) and transverse (r±) direction with respect to the ep cms, our chosen 
global frame of description. The left (right) panel corresponds to small-x (large-x) events 
with fixed Q 2 = 28 GeV 2 . In each plot the three curves correspond to times 0.4/12/20 
fm/c after the photon-quark scattering, and each curve includes all particles (partons and 
formed hadrons), which are actively present in the mixed particle system at the specified 
times. Comparing the time development of the spectra dN/dy spectra for small-x and 
large-x events, one observes that in the former event type most particles are produced at 
central rapidities \y\ < 1, with a shift toward the proton side (negative y), while in the 
latter event class this region is least populated. Related to that, the dN/dz distributions 
show a much larger particle production along the beam axis around z = fm for small-x 
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than for large- x events. ^From the l/NdN/dr\ distribution one can see that for small- 2; 
scattering the diffusion of the expanding particle system in the transverse direction is faster 
than in the large- x events, an effect that arises from the transverse pressure of the larger 
number of produced particles in the central region. 

4.2 Inclusive hadron spectra in xp, p^, and the (p^ 2 ) dependence 

The study of particle multiplicities and momentum distributions of the hadronic final 
state at HERA provides sensitive information about both the QCD processes at the parton 
level and the properties of hadron formation. An excellent recent review can be found in 



Ref. [41|. One of the attractive features of the HERA experiments is the production of a 
large-mass hadronic final state with W — 100-250 GeV (an order of magnitude larger than 
in previous fixed-target experiments). The conjecture is therefore that the influence of the 
non-perturbative QCD effects is less important and that - in the spirit of 'local parton- 
hadron duality' |12| - the observed hadronic final state reflects more the dynamics of the 
partonic processes. 

Particularly sensitive measures of the parton level dynamics are the x p and p\ distri- 
butions, as well as the (pj_ 2 ) of produced charged hadrons, as measured in the jp cms (||), 
where the Feynman variable xp = 2pt/W and p*±_ characterize the momentum components 
of hadrons parallel and transverse to the photon direction. At large values of W (small 
values of Bjorken x), these observables are sensitive to hard multigluon radiation. This 
feature is evident in Fig. 8, where we plot our model results for the xp and p\_ spectra and 
the dependence of (p*i 2 ), for three typical HERA values Q 2 = 8/14/28 GeV 2 and the small- 
x regime defined by (p7|). Also shown are the corresponding measured distributions for 
Q 2 = 28 GeV 2 from ZEUS @, with which the calculated dashed-dotted curves (Q 2 = 28 



GeV 2 ) agree reasonably well. All three distributions have a specific form due to QCD gluon 
emission on the parton level which cannot be explained by the naive 'quark parton model' 
(QPM) which accounts only for the lowest-order photon-quark scattering and omits all 
higher-order QCD radiation. For comparison, the QPM results are plotted as thin curves. 

The xp distributions (top) show steep exponential decreases above xp ~ 0.05 — 0.1 
and an enhanced particle yield below that value. Note that the QPM, i.e., the leading- 
order Born scattering alone, gives a slightly shallower decrease. The effect of the higher- 
order radiative processes is however very prominent in the p*±_ spectra (middle) integrated 
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over xp > 0.05, which show a power-law dependence due to multi-gluon emission and a 
significant contribution of hard gluons with transverse momenta > 3 GeV (large, in view of 
vQ^ = 3-5 GeV). This result is in vivid contrast with the corresponding QPM result, which 
hardly gives any transverse momenta ^ 1 GeV. The mean square of p\, {p*i_ 2 ) (bottom) is 
particularly sensitive to the tail of the distribution and exhibits half of a 'seagull' shape for 
positive values of xf- The rise of (p± 2 ) with increasing xf is due to the leading hadron 
effect, which can be understood as follows: If denotes the fraction of momentum the 
initially-struck quark transferrred to a hadron, 

(P*± 2 ) = *h ( (Pl 2 )prim + {pl^seej + (p\ 2 ) frag , (38) 

where p*± pr i m is the primordial transverse momentum of the quark to which the photon 
couples (c.f., eq. (|?])), p± S ec is the secondary contribution from QCD radiation, &ndp± f rag 
denotes the additional transverse momentum produced in the fragmentation and hadron 
formation process, which is in the average about 0.45 GeV and almost independent of W and 
Q 2 . /,From ( p8| ) one sees that the leading hadrons at large xf, which carry higher fractional 
momentum z^, also carry a higher fraction of the actual parton transverse momentum 
(primary plus secondary). This explains qualitatively the 'seagull' shape and the rise of the 
(p*L 2 ) of hadrons as a function of xp. Furthermore, one observes that, the more gluons with 
relatively large transverse momentum are radiated, the larger the contribution to (pj_ 2 }secj 
and hence the stronger is the effect. In the QPM with no QCD radiation at all, the total 
(p*L 2 ) is therefore a factor of 5-10 smaller and has only a weak xf dependence. 

In Fig. 9 the W and Q 2 dependences of the mean squared p\_ of hadrons is shown for 
two intervals 0.1 < xp < 0.2 and 0.2 < xf < 0.4, and compared with data obtained by the 
ZEUS collaboration. The agreement with the data is fairly good for the <3 2 -dependence, 
whereas it is less clear for the dependence on W. The (pj_ 2 ) depends strongly on both 
W (for fixed Q 2 = 28 GeV 2 ) and Q 2 (for fixed W = 120 GeV). It is worth noting that 
in previous fixed-target experiments at lower energies the (p*^ 2 ) is generally much smaller, 
depending only weakly on W (however at much smaller W values), and essentially flat in 
Q 2 - 

We remark that the new class of diffractive events with a large rapidity gap is measured 



[ 43 1 to have very different p*j_ spectrum and ) from the non-diffractive events that we 
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have just discussed. The particle distribution in p* ± of diffractive events falls much steeper 



are smaller by a factor of 2-5, and lie just slightly above the QPM curve in the bottom plot 
of Fig. 8. This indicates that diffractive events resemble in a way DIS events with very 
little QCD radiation, consistent with the common interpretation that the photon couples 
in these events to the proton via a colourless intermediate state which does not fragment 
by multiparton emision. 

4.3 Transverse energy flow 

Whereas the inclusive xf and p± distributions of produced hadrons are sensitive to the 
multi-jet structure due to hard-gluon radiation as discussed above, the analysis of inclusive 
hadron distributions in terms of the global energy flow extends to classes of events which 
cannot be identified unambigously as n-jet events. According to the idea of 'local parton- 
hadron duality' |^], the pattern of overall distribution of energy among the partons in 
an event determines the energy flow observed in the hadronic spectrum. The energy flow 
dE^/dr]*, and similarly the particle flow dNh/drj*, are commonly studied as a function of 
pseudorapidity rf in the hadronic centre-of-mass system, i.e., the "fp cms, where 



In Fig. 10 we show model results for the distribution of hadrons in rj*, E*^_ and as 

well as the E^_ flow. As before, we choose Q 2 = 8/14/28 GeV 2 for the small- 2; regime 

2.3 • 1(T 4 < x < 1.7 • 1(T 3 corresponding to W > 60/90/130 GeV (c.f. Table 1). The 

plotted distributions reflect the typical event geography that one expects already from QPM 

considerations. Recall that in the QPM, with the neglect of higher-order QCD corrections, 

the struck quark and the proton remnant system each carry an energy W/2 in the "yp 

cms and move back to back with rapidities ±yj^ ai . oc ± In W/m^ . The fragmentation 

of the two receding charges fills the intermediate pseudorapidity region with hadrons]^]. 

The width of the hadron distribution in the final state is proportional to InVF, while its 

height is approximately independent of W. The width of the quark jet and the proton 
10 For the purposes of this discussion, we will now neglect the difference between rapidity y and pseudo- 
rapidity rj. 



and resembles closely the QPM curve in the middle plot of Fig. 8. The mean values (p^_ 2 ) 





(39) 



with 9*,p*^,p* L defined with respect to the photon direction. 
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fragmentation region bounded by -±.rf max is typically about 2 units, corresponding to the 
x^-range \xp\ > 0.05. Hence, at high values of W (small x), the pseudorapidity range 
populated by hadrons can be divided in three regions: (i) the current jet region from 
{ r lmax ~ 2) to rf max , (ii) the proton fragmentation region from [—rftnax + 2) to —rf max , and 
(hi) a central plateau region in between. 

The pseudorapidity distribution dNh/dn* (Fig. 10, top), as calculated in our model, 
shows a distorted version of the naive QPM picture due to the higher-order QCD radiation 
effects. The spectrum is asymmetric and the central plateau is rising from the proton 
fragmentation region to the current jet region, rather than being flat. Particularly different 
from the QPM is the behaviour in the current jet region rf ^ 3, which shows a clear increase 
with Q 2 of both the height and the width of this part of the hadron distribution, an effect 
of the jet broadening due to gluon radiation off the struck quark. 

The transverse energy distribution dNh/dE*^ (Fig. 10, middle) is naturally similar to 
the p* ± spectrum discussed before (c.f. Fig. 8). It again exhibits a power-law behaviour 
that is characteristic for gluon emission with significant p*^, a feature which becomes more 
prominent with increasing Q 2 , because of the enlarged phase-space and extended duration 
of parton shower activity before hadronization. 

The hadronic energy flow dE^/drf (Fig. 10, bottom) mirrors the distribution of energy 
and transverse momentum among the final-state particles in a similar way to the pseudora- 
pidity distribution discussed above. The characteristic features are: first, a central plateau 
with a slight dip and a height almost independent of Q 2 , secondly, an increase with Q 2 
of energy deposit in the current jet region around rf ~ 3, resulting from radiation of the 
time-like shower of the quark after the hard scattering, and thirdly, a similar though much 
less significant increase with Q 2 of the activity around rf ~ —5 in the proton fragmenta- 
tion region due to radiation from the space-like shower before the hard interaction, which 
is going along the proton direction and opposite to the time- like radiation. 

In Fig. 11 the x (or, equivalently, W) dependence of the Ej_ flow in the current jet 
hemisphere in the ^yp cms is investigated for fixed Q 2 = 28 GeV. We plot dE^/drf for three 
different x ranges with average values (x) = 3.7- 10~ 3 /1.5- 10 _3 /7.2- 10 -4 , corresponding to 
W ~ 78/137/197 GeV. Also depicted are the data points of the measured distributions from 



the ZEUS collaboration [44]. Although the model slightly underestimates the data around 
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rf = 0, qualitative conclusions that may be drawn are: first, the height of the plateau-like 
region for r/* < 1-5 is rather independent of x, and secondly, with decreasing (x) (from top 
to bottom) the peak around the current jet moves visibly towards larger rapidities, with 
Vpeak — 2.4/2.8/3.2, respectively, while the height of the peak appears to be stable. 

4.4 Mass distributions of the observed hadronic final state 

As mentioned in the introduction, a new class of DIS events is observed at HERA 
in the experiments by ZEUS |B[] and HI ||, events which are characterized by a large 
rapidity gap (LRG) between the proton and the rest of the hadronic final state that is 
measured in the detector. The properties of these events indicate a diffractive production 
mechanism via exchange of a coherent colourless object between the photon and the proton 
(c.f. Fig. lb), accompanied by a suppression of QCD radiative processes, which are, as 
we discussed, so prominent in non-diffractive events (c.f, Fig. la) with no rapidity gap 
(NRG). Because in the experiment both diffractive LRG and non-diffractive NRG events 
are mixed (with a relative contribution of ~ 5-10 % from LRG events), the determination 
of the diffractive cross-section requires the detailed knowledge and subtraction of the non- 
diffractive contribution. 

A method to separate diffractive and non-diffractive contributions suggested by ZEUS 
[fj5f uses the mass Mx of the hadronic system X that is measured in the detector, where 

M x - | " (ix) - (ix) - (ix) l < W\ (40) 

includes all observed particles except for the outgoing electron. Because of the finite res- 
olution and geometric acceptance of the detector, Mx is naturally smaller than the total 
invariant mass W, with the event fluctuations giving rise to a distribution in Mx- The 
remarkable feature of the distributions in M x and In M x is that they exhibit very different 
behaviour for the two event types and are sensitive measures of the event structures. In this 
context, we investigate in the following the Mx spectrum of purely non-diffractive NRG 
events and its Q 2 and x (or W) dependences, so to provide an estimate of the non-diffractive 
contribution underlying the diffractive LRG component. 

Let us briefly summarize the state of knowledge in order to set the stage. As illustrated 
in Fig. lb, in diffractive scattering the outgoing proton or low-mass nucleonic system 
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remains colourless and escapes through the forward beam hole, while the system X from the 
dissociation of the photon is, in general, almost fully contained in the detector. Diffractive 
dissociation prefers small Mx values (In M x <, 4) with an event distribution of the form 

© 

JKfdiff ( 1 \ n 

(41) 



dhxM x " \M X 



where at high energy one expects n ~ [ 40 1 , approximately independent of the total 
cm energy W. 

On the other hand, in non- diffractive events the incident proton is broken up and the 
remnant of the proton is a coloured object with the struck quark taking away the net 
colour. As we discussed before, this results in a substantial amount of initial- and final- 
state radiation followed by hadron formation between the directions of the proton and the 
current jet (c.f., Fig. la). ^From perturbative QCD arguments, as well as simple phase- 



space considerations, one expects [45| that the associated event distributions are peaked at 
large Mx values (lnMj^ ~ 5-10) with an exponential fall-off towards smaller masses, P] 

( jj l j-nondiff 



d\nM x 



c exp (b lnM|) , (42) 



where c is a constant. The slope b is the parameter of interest: on the parton level it can 
be shown to be determined by the QCD Sudakov form factor and thus the probability for 
gluon emission. In our model for parton-hadron conversion, it is therefore closely related 
to the probability for cluster formation and hadron production. 

In Fig. 12 we show the non-diffractive event distributions in Mx (top) and lnM^ 
(bottom). The plots compare calculations with the fixed values Q 2 = 8/14/28 GeV 2 , 
normalized to the total number of events, and within the nominal kinematic acceptance of 
the ZEUS detector, -3.8 < r) < 4.3. The distribution 1/MdM /dM x (M = M nondiff ) in the 
top part of the figure exhibits a clear Q 2 dependence in both the position of the peak and 
the extension of the tail towards large Mx values. The mean values (Mx) = 78/93/105 
GeV for Q 2 = 8/14/28 GeV 2 . The properties of the distributions are most evident when 
studied as a function of lnM^, as in the bottom part where 1/NdN /d\w.M x is shown for 
11 Another salient feature of the lnM^ distribution is an observed scaling in \jxMx — lnV^ 2 jlB]], im- 
plying that the position of the high-mass peak in InMx grows proportional to lnW^ 2 , and the slope of the 
exponential fall-off to small lnM^ values is approximately independent of W. 
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the same parameters. In this representation the the mass peak exhibits a steep exponential 
fall-off towards smaller M x values. The associated slope exhibits no significant dependence 
on Q 2 and comes out as b = 0.95 ± 0.1, when compared with the form ( f42|) as indicated by 
the straight line in the plot. On the other hand, an experimental data analysis by ZEUS 
|45f yields a steeper slope, namely b exp = 1.46 ± 0.15. The discrepancy between b exp and 
b can have various reasons associated with experimental effects that we did not attempt 
to simulate or correct for, e.g., detector acceptance or other effects, such as energy loss of 
particles in the calorimeter. Such effects may affect the value of b. However, the fact that 
both b and b exp come out to be universally constant supports the the conjecture that the 
difference between the values of b and b exp is due to global effects that are missing in our 
calculations. 

In Fig. 13 we investigate the W dependence of the In M x distribution at Q 2 = 14 GeV, 
plotting our results for three distinct intervals of the total 7p cm energy W. One observes 
that the slope is the same in all three W ranges, and hence appears to be independent 
of W as well as Q 2 . The position of the peak, however, is shifted to larger values as W 
is increased, as in the previous figure when Q 2 is increased. Also shown in Fig. 13 are 



the corresponding measured distributions measured by ZEUS [45], with which our model 
calculations agree reasonably well for the large Mx range, lnMj^ > 4 — 6, but which show 
an additional component at small values In M x < 4. This latter is due to the diffractive 
(LRG) contribution, which evidently has a plateau-like (rather than exponential) shape, 
in agreement with the expectation (fil|). The comparison between the data points and our 
model results exhibits two important conclusions. First, the diffractive component cannot 
at all be explained by the standard QCD parton shower evolution plus hadronzation model 
of non-diffractive events: in fact, it is completely absent therein. Secondly, the diffractive 
and non-diffractive contributions appear to be sharply separated when studied with respect 
to the variable ]xlM x , and allow one to subtract cleanly from the measured data sample 
the non-diffractive part, as calculated using this or other QCD parton shower models. 

Finally, we would like to comment on the difference between the value of b calculated 
within our model, as compared to other parton shower models [19, 21, 22 1, which use 
the string fragmentation approach [Eq] to hadron formation from the final-state parton 



ensemble. As investigated in Ref. [45|, the latter give a value of b ~ 2, i.e. about twice as 
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large as in our model. We believe that this difference arises from differences in modelling 
the parton-hadron transition, i.e., string fragmentation versus cluster formation and decay 
PI , as we now discuss in more detail. 

Let us first recall some general features of hadron distributions and correlations within 
the Mueller approach [fj7f| , where they are related by unitarity to appropriate absorptive 
parts of forward multiparticle scattering amplitudes. In the beam fragmentation region, 
which is relevant to this discussion, asymptotic properties of the single- and multiparticle 
distributions are controlled by Regge singularities. In particular, the asymptotic value of 
the single-particle density is controlled by the pomeron, with subasymptotic corrections and 
finite-range multiparticle correlations controlled by subleading Regge singularities. These 
give in particular local two-particle correlations with a correlation length Ay = 1/Aa, 
where Aa is the difference between the intercepts at t = of the pomeron and the next 
subleading trajectories, commonly believed to be the p and degenerate trajectories with 
Aa ~ 1/2, yielding Ay ~ 2. The it trajectory with Aa ~ 1 would yield shorter-range 
correlations with Ay ~ 1. These subleading trajectories would also yield subasymptotic 
corrections to the Mx distribution: b ~ 2Aa, corresponding to b ~ 1(2) for the p(tt) 
trajectories, with the tail corresponding to rapidity-gap events corresponding to pomeron 
exchange with b ~ 0. 

With these points in mind, we now recall aspects of particle production according to 



the Lund string fragmentation model [46], as compared to our approach. In the former 
model, the string is chopped at an ordered sequence of points along the rapidity axis, with 
separations chosen randomly but with a mean value 6y ~ 1. The string bits then hadronize 
independently, with resonance decays resulting in a correlation length Ay ~ Sy ~ 1. The 
adjacency in rapidity of the Lund fragmentation model clearly results in the minimum 
possible correlation length Ay, and hence effectively to the largest possible Aa. Thus it 
is no surprise to find that simulations based on this model yield a relatively high value of 
b ~ 2, corresponding to the tt trajectory in the Mueller language. On the other hand, in our 
space-time approach, the pre-hadronic clusters are formed by adjacent pairs of partons in 
position space, which are not necessarily the closest in rapidity space. This point is reflected 



12 The preceding parton shower stage is essentially the same in our and other models ]lq| , with our 
additional space-time information becoming relevant during cluster formation and hadronization. 
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in Fig. 14, where we see that the separation in rapidity between partons that combine to 
form a cluster is typically 5y ~ 2, approximately a factor 2 larger than in the Lund model. 
We therefore expect in our model that Ay ~ 2, corresponding to Aa ~ 1/2 and b ~ 1 as 
we found above. 

To the extent that the experimental value of b exceeds unity, it may be that our space- 
time approach deviates too far from the adjacency in rapidity of the Lund string fragmenta- 
tion model, and the truth may well lie somewhere in between, corresponding in the Mueller 
approach perhaps to a combination of the p and n trajectories. One way to test this would 
be to measure experimentally the rapidity correlation length, and compare it directly with 
the predictions of various models. An interesting issue to watch will be whether b and the 
effective two-particle correlation length depend on Q 2 or W 2 . The naive Mueller Regge 
described above has been derived for incident hadrons, and may require modification at 
large Q 2 . 

5. SUMMARY 

We have presented in this paper the application to DIS at HERA of a model for the 
quantum kinetics of multiparticle production that includes the space-time development 
of the parton shower, cluster formation and hadronization. Compared with our previous 
work, novel features include tracking back to the initial proton the development of the 
space-like parton shower prior to its interaction with the virtual photon radiated by the 
electron. Our procedure tracks in space and time the emission and evolution during this 
development of time-like secondary partons, as well as the spectator partons in the proton 
beam fragment]^]. As in our previous work, the coalescence of partons to form pre-hadronic 
clusters is determined statistically by a spatial criterion motivated by confinement and a 
simple non-perturbative model for hadronization. 

Our space-time approach has enabled us to map the history of the particle densities 

and associated spectra, including the rapidity, longitudinal and transverse distributions of 

particles. These hard results may be compared with intuitive pictures of the space-time 
13 This machinery will be applied in the future to more complicated situations including eA, pp, pA and 
AA collisions. 
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development of hadronic final states in DIS. They will also form the basis for the subsequent 
extension of our approach to shadowing and other interesting effects in eA scattering. 

We have also explored in our model inclusive hadron spectra in xf and px, and the 
transverse energy flow. Although our model reproduces the general features of the observed 
pattern in energy flow, it shares with other simulations the tendency to undershoot the data 
around 77* = 0. However, the discrepancy is not dramatic, and does not make a strong case 
for the presence of imortant physical effects not present in the MLLA approach we use here. 

Our model provides distributions of the mass Mx of the observed hadronic final state 
in events without a large rapidity gap, which can be used to estimate the background to the 
cross-section for LRG events. We find a spectrum ~ exp(blnM|-) with an exponent b ~ 1, 
which is not very different from the observed value 6 exp ~ 1.5. A detailed comparison with 
the data requires more understanding of detector effects and final-state hadron interactions, 
which goes beyond the scope of this paper. The value of b is sensitive to the rapidity density 
and other properties of pre-hadronic clusters, so the relative success of our model, which 
has no parameters adjusted from its previous applications to e + e~ annihilation, gives us 
some hope that we are capturing important aspects of this physics. 

As already mentioned, the new features of our approach introduced in this paper, in- 
cluding the space-time treatment of the initial hadronic state, open the way to future 
applications of our model to eA, pp, pA and AA collisions, where the novel features associ- 
ated with high parton densities will become more marked. We aim eventually at a unified 
space-time description of all these hadronic processes. 
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TABLE CAPTIONS 



Table 1: 

Examples of the kinematic relations (g), (g) between the Bjorken scaling variable x, the 
the absolute squared invariant mass of the photon Q 2 , and the inelasticity variable y, as 
well as the total invariant mass W of the hadronic system. 

FIGURE CAPTIONS 
Figure 1: 

Schematic diagram of particle production in a) non-diffractive, and b) diffractive DIS 
events. Here W is the total invariant mass of the produced hadronic system, Mx is the mass 
of the observed final state in the detector, and in b) P' represents the outgoing (excited) 
proton or low-mass nucleonic resonant state. 

Figure 2: Deep-inelastic ep scattering viewed in a) the HERA laboratory frame where e 
and p collide head on, and b) the 7p cms, where 7 and p collide head on. 

Figure 3: 

Contour plot in the x-Q 2 plane of the mass of the produced hadronic system W and 
the inelasticity variable y for y/s = 296 GeV at HERA. The approximate region of previous 
fixed-target experiments is indicated by the shaded area at£>l(T 2 . (FromRef. El.) 

Figure 4: 

Schematics of the components of our model for DIS: The highly Lorentz-contracted 
incoming proton with its initial-state parton configuration evolves from the remote past 
t = t < 0, and is struck by the photon at t = 0. This hard interaction picks a quark 
out of the proton's parton cloud, thereby triggering initial-state (space-like) and final- 
state (time-like) parton showers. With increasing time t — * +00, the partons evolve by 
further radiation, whereas the remnant proton propagates on as a coherent remainder. 
In the process of hadronization the produced partons may coalesce to colourless clusters 
if they are nearest neighbours in space-time, whereas the virtual partons of the proton 
remnant combine with a colour- neutralizing parton to form a massive beam cluster. Both 
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'parton clusters' and the 'beam cluster' subsequently convert into primary hadrons that 
subsequently decay to low-mass final-state particles. 



Figure 5: Schematic diagram of the probabilistic parton evolution in the MLLA framework 
(solid lines are quarks, curly lines are gluons). The initial-state quark with space-like 
virtuality — p\ ~ Qq evolves from t = to forward in time and toward the hard ^yq vertex by 
successively increasing its off-shellness up to —p^ ~ Q 2 , when it is struck by the photon 
at t = 0. The outgoing quark is provided by the momentum transfer with a time-like 
virtuality ~ Q 2 , and radiates off its excitation by successive gluon radiation until it starts 
to hadronize by coalescence with another parton, at which point the shower terminates 
naturally. 

Figure 6: 

Characteristic differences in global properties between 'smail'-x (left panel) and 'large'- 



x DIS events (right panel), as defined in (37). The two ranges refer to the values of Bjorken 
x of the quark struck by the photon with selected Q 2 . The small- x range is typical for the 
kinematics of HERA experiments, whereas the large- x regime corresponds to the phase- 
space region probed by previous fixed-target experiments. Compared are a) the probability 
distributions for the production of a hadronic system of mass W (top), b) the corresponding 
W dependence of the total hadron multiplicity Nh (middle), and c) the resulting mean 
values of xp = 2p z /W, with p z in the ep cms along the beam axis in the opposite direction 
to the incoming proton. 

Figure 7: 

Characteristic differences in the space-time evolution pattern of 'small'-x (left panel) 
and 'large'-x DIS events for DIS events at Q 2 = 28 GeV (in correspondence with Fig. 6). 
Compared at 3 different times t = 0.4/12/20 fm/c are a) the rapidity distribution dN/dy 
(top), b) the particle distribution along the ep beam axis, dN/dz, and c) the particle 
distribution perpendicular to the beam axis, l/NdN/dr±. Note that the distributions 
include all particle species (partons and hadrons) present in the system at the quoted given 
times. However, the y distributions only count the secondary particles and exclude the 
primary partons of the original proton, beacuse their rapidity is not well defined, whereas 
the z and r± spectra include also those primary partons. 
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Figure 8: 

Model results for differential charged hadron multiplicities with respect to the 'yp cms 
as a function of xf and p*^_ for xf > 0.05, as well as {p*^ 2 ) as a function of xf- We compare 
results for the small-x regime (g7|) at typical HERA values of Q 2 = 8/14/28 GeV 2 . The 
data points are measured distributions from ZEUS |^3| for Q 2 = 28 GeV 2 , and the thin 
solid curves represent the corresponding expectations of the 'naive' quark parton model 
(QPM). 

Figure 9: 

Model results for the W dependence for fixed Q 2 = 28 GeV 2 and the Q 2 dependence 
for fixed W = 120 GeV of the mean squared transverse momentum of charged hadrons 
{p*i_ 2 )- The plots refer to the 'yp cms and separate the two intervals 0.1 < xp < 0.2 and 



0.2 < xf < 0.4. The data points are from the ZEUS experiment 43] 



Figure 10: 

Hadron distributions in the "fp cms as a function of pseudorapidity rf and of transverse 
energy Ej_, as well as the E*^_ flow versus pseudorapidity. The model results refer to Q 2 = 
8/14/28 GeV 2 and 2.3 • 10~ 4 < x < 1.7 • 10" 3 corresponding to W > 60/90/130 GeV. 

Figure 11: 

The x (W) dependence of the E^_ flow on the current jet side in the jp cms at Q 2 = 28 
GeV and for (x) = 3.7 • 10~ 3 /1.5 • 10" 3 /7.2 • 10~ 4 , corresponding to W ~ 78/137/197 GeV. 
The experimental data are from ZEUS 



Figure 12: 

Comparison of normalized event distributions in Mx (top) and In M x (bottom) at Q 2 = 
8/14/28 GeV 2 , and for 2.3 • 10~ 4 < x < 1.7 • 10~ 3 . with -3.8 < i) < 4.3. 

Figure 13: 

Model results for the W dependence of the In M\ distribution at Q 2 = 14 GeV for three 
W intervals, compared with the corresponding measured distributions from ZEUS J45jJ . 

Figure 14: Top: Population of rapidity y* in the jp cms of pre-hadronic clusters formed 
from coales — ced parton pairs in the current jet region. Bottom: Distribution in relative 



41 



rapidity 5y* = y* — y\ of parton pairs making up the clusters. The thin full line represents 
a constant rapidity separation Sy* = 1 and serves as reference to the discussion in the text. 
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